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SUMMARY 

A thermo-elastic  model  of  a Laser  irradiated  plate  is  developed. 
The  model  accounts  for  plate  flexure  and  shear  deformation  as  well  as 
plate  stretching.  An  exact  solution  in  series  form  is  found  for  the 
dynamic  response  of  a circular  plate  clamped  at  its  boundary.  The 
plate  is  subjected  to  Laser  irradiation  at  its  center,  and  normal  to 
the  plate  surfaces.  Experiments  with  a Holobeam  model  630-Q  Nd  Laser 
were  conducted  to  verify  the  mathematical  model.  A careful  comparison 
between  theory  and  experiment  of  transverse  plate  deflection  induced 
by  the  Laser  shows  good  agreement. 
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I 

I 

1.0  INTRODUCTION. 

It  is  well  known  that  Lasers  have  the  Inherent  ability  to  deposit 
focused,  radiant  energy  on  a structure  or  its  elements.  Lasers  can 
concentrate  a short  duration,  high  energy  flux  in  an  extremely 
narrow  beam.  The  effect  of  short  duration,  high  energy  irradiation 
on  an  opaque  solid  can  take  several  forms  which  are  described  (roughly) 
as  follows  (see  Chapter  3 of  ref.  5) : 

(a)  Conplete  local  vaporization  of  the  material,  and  the 
resulting  creation  of  openings  (holes).  If,  in  addition, 
the  structure  so  punctured  is  in  a state  of  initial  stress, 
there  will  be  an  additional  dyneunic  effect  due  to  unloading 
waves  and  wave  reflections,  with  the  possibility  of  stress 
intensi f ication . 

(b)  Sudden  deposition  of  thermal  energy,  without  a change  in 
phase.  This  causes  sudden  thermal  stresses  in  the 
structure,  and  because  of  the  rapidity  of  the  energy 
deposition  process,  there  will  be  thermally  generated 
stress  waves. 

(c)  It  is  possible  for  the  structure  to  experience  partial 
surface  vaporization  over  the  effective  Laser  beam  cross- 
section.  This  results  in  material  removal  and  plasma 
generation.  Subsequent  heating  of  the  plasma  gives  rise 
to  shockwave  formation,  resulting  in  impulses  transmitted 
to  the  solid. 


Effects  (a),  (b),  and  (c)  can  coexist,  but  they  can  be 
separated  for  analytical  purposes.  The  present  analytical  and  experimental 


4 


investigation  is  concerned  with  case  (b) , above.  We  consider  the 
specific  case  of  a thin,  circular  plate,  rigidly  clanqped  at  its 
boundary.  The  plate  is  subjected  to  Laser  irradiation  at  its  center, 
normal  to  the  plate  surface.  The  beam  energy  density  is  adjusted  to 
cause  heating  of  the  plate,  but  to  avoid  phase  changes  (melting)  of 
the  plate  material.  Our  objective  is  (a)  to  develop  a mathematical 
model  which  predicts  the  thermoelastic  dyn2unic  response  of  the  plate, 
and  (b) , to  perform  suit2d3le  experiments  to  establish  the  validity  of 
the  mathematical  model. 
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2.0  THERMO-ELASTIC  PLATE  THEORY. 


2.1  Plate  Strain  and  Stress  Components. 


We  shall  derive  a themoelastic  plate  theory  which  accounts  for 
stretching  and  bending  of  the  plate.  With  regard  to  stretching,  it 
is  assumed  that  the  plate  stretches  uniformly  through  the  const2uit 
plate  thickness  h,  with  components  u^Cx^jX^.t).  With  regard  to  plate 
bending,  it  is  assumed  that  straight  material  elements  (lines) 
originally  normal  to  the  plate  median  surface  remain  straight  but 

ate  through  an  angle  following,  Greek  subscripts 

ive  on  the  range  1,2;  while  Latin  subscripts  assume  values  1,2,3. 

Thus  the  displacement  vector  u^,  i=l,2,3  assumes  the  form 


»XT»t)  + u (x  ,x  ,t) 
al2  a 12  al2 


u (x, ,x,,z,t)  *w(x, ,x_,t) 
Z i.  2 1 Z 


The  general  strain-displacement  relation  is 


E ••  — (u  + u ) 
‘^ij  2 '"i,j  '*j,i' 


and  upon  substitution  of  (1)  into  (2)  we  obtain 


'.6  ■ T “6,«’ 


e - e - —('I'  +w  ) ; e ? o 
oz  za  2 a ,a  zz 


It  will  be  convenient  to  define  the  plate  strain  components 


12  h/2  12  h/2 

-h/2  h^  -h/2  ^ ^ 
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, h/2 

^ ^ -h/2 


h/2  h/2 

f e dz  = — / e dz 

■h/2  ^ -h/2 


1 h/2  ^ h/2 


similarly,  we  define  the  plate  stress  components  by 


h/2  h/2 

/ T oZdz  = / T„  zdz  * M, 

J raO.  * Qrt  h 


h/2  h/2 

/ T dz  = J 7 dz 
cxz  ^ zot 

-h/2  -h/2 

h/2  h/2 


In  view  of  (3)  and  (4),  the  plate  strain-displacement  relations  are 


“aP  2^'^a,P  * 


a * it  -*•  w 
^ a ,a 


"as  ■ I'Ve  * -e.a’ 


and  we  note  that 


e.-zma-J-n-  , £ ■£  “-rq^ 

aP  aP  aP  az  za  2 ^ 


2.2  Stress  Equations  of  Motion. 

Me  shall  now  derive  the  pertinent  equations  of  motion  with  the 
aid  of  Hasillton's  principle 


t 
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(6k  - 6u  + 6w)dt  - 0 
t. 


(8) 


where  6k  and  6u  are  the  first  variations  of  kinetic  and  potential 
energy  in  the  plate,  respectively,  and  6w  is  the  virtual  work  of 
external  forces  in  an  (admissible)  variation  of  the  plate  displacement 
field.  The  kinetic  energy  of  a solid  with  volume  V is  given  by 


K “ -J  / J / p(Uj^  “2 


V — V 

With  reference  to  (1)  we  have 


(9) 


II  = zw  + u 
^ a 


u ^ w 
3 


(10) 


2md  upon  substitution  of  (10)  into  (9) , we  obtain 
1 2’  * 2 

K = j / J P(2  + w + u^u^+  2zu^^i^^)dz  dA 

A ~h/2 

- j / (phw^  + — + phu^u^)dA  (11) 

A 

Taking  the  first  variation  of  (11)  with  respect  to  displacements, 

6k  « / (phw6w  + ph\  + phu  6u  )dA  (12) 

i 12  ^ a a a 

A 

We  now  integrate  both  sides  of  (12)  with  respect  to  time  from  t^  to  t^, 

and  then  use  integration  by  parts  on  the  right  hand  side.  In  line 

with  the  requirements  of  Hamilton's  principle,  we  set  6ij»^**6w-6u^«0  at 

t>t,  and  at  t>t„.  Thus  we  obtain  the  result 
1 2 


6xdt  ■ ” /^  / (phw6w  + * phu^6u^)dAdt 

ti  tj_  A 


(13) 
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The  variation  of  the  strain  energy  density  of  the  plate  is  given 
by  (see  (37)) 


6^  = M „6in  o + Q 6q  + N -6n  „ 
P oB  aB  a a oB  aB 


Consequently,  the  variation  of  the  total  strain  energy  is  given  by 

6U  = / S'hr  dA  *=  / (M  dm  o + Q 6q  + N „6n  „)dA 

^ p ^ aB  aB  a ^a  aB  aB 

= / (M  o+Q<5iJ»  + q6w  +n  _6u  -)dA  (IS 

^ aB  a,B  a a a ,a  aB  a,B 

where  we  have  used  (6) • Upon  application  of  Green ' s theorem  ir.  the 
plane  to  (15) , we  obtain 

So  . |<M;,gS<,^i.g  . i^Swi,^  * 

-[  ' '"aB.B-^a’^^a  " * '‘aB.B®“a"“ 

where  C is  the  boundary  cur^'e  of  the  simply  connected  domain  of  the 

A 

plate,  and  n,£  are  unit  vectors  which  are  normal  and  tangential, 
respectively,  to  C.  We  also  note  the  relationships 


u 

m 

n u 

Ub  • 

£ u 

1 

n 

a a 

£ 

a a 

'^n 

m 

"a*a 

£ 

a’^a 

On 

m 

"A 

0£- 

Va 

» 

M 

nn 

m 

n n-M  a • 

a B aB 

“n£- 

Vb^oB 

N 

nn 

m 

n n„N  a 

a B aB 

•*n£- 

Vb^oB 

"aB*VB  ■ "nn**n  * 

N -6u  n-  ■ N 6u  ♦ H ,6u, 
aB  a B nn  n n£  t 
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It  can  be  shown  (p.  465,  eq.  (14.19),  ref.  [1])  that  the  virtual 
viork  of  external  applied  forces  and  boundary  forces  is  given  by 


6w  - / p6wdA  + f (M  M 

^ ^ i\n  n nJC  X*  n lui  n nx»  JC 


(19) 


where  p is  the  transverse,  normal  pressure  acting  on  the  plate  in  A 

* * 

and  M , Q , etc.,  are  stress  resultants  acting  on  the  boundary  C. 
nn  n 

Upon  substitution  of  (13) , (16) , and  (19)  into  (8)  we  obtain 

f / 1 <-PhiJ+Q„^„+P)  6w+  (-  ^ g-Qa’  % 

+ (-PhVN^e^g)6u^ldA 

. /2 


+ (N*.-N  -)6u.  + (n*  -n  )6u  Idlldt  - 0 
n£  nr  £ nn  nn  n 


(20) 


Since  A,  C and  ^2~^1  arbitrary,  the  two  integrals  in  (20)  vanish 
separately.  Furthermore,  the  displacement  variations  are  arbitrary 
in  A and  on  that  portion  of  C where  stress  resultants  are  prescrDsed. 
We  thus  infer  that  in  A 


°a,a''P  “ ° 
" U ”a6,e 

^ ”a6,B  ’ ° 


Q„-0 


(21a) 

(21b) 

(21c) 


emd  on  C one  member  of  each  of  the  following  pairs  is  specified: 


(w. 

V- 

"nn 

<“»■  "nt 

»*n£> 


(22a) 


(22b) 
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2.3  Constitutive  Relations  and  Displacement  Equations  of  Motion. 

We  assume  that  the  Neumann-Duhamel  relation  expresses  the  relation 
between  stress,  strain,  and  temperature  in  a homogeneous,  isotropic 
material,  i.e. , 


EE.  . * (l+v)T,  . - VT  + EaT6.  . 

ID  Id  ij  ID 


For  our  present  purpose  we  wite 


EE  o = (1+V)T  o-VT  6_-vr  6„  + EtxTfi  „ 

ae  ae  YY  zz  oB  aB 

T « 2(C^Ge 
za  za 


and  we  ignore  the  relation 


(23) 


(24a) 

(24b) 


EE  - T - VT 
zz  zz  YY 


If  we  set  o-^B  in  (24a) , 


E + 


oa  (1-v)  oa 


and  therefore 


+ E<XT  * 0 

w obtain 

2v 

(1-V)  ^zz 


2EaT 

(1-v) 


(24c> 


(25) 


T m - . vr»  r A V X CAJLX  ^ 

aB  “ (l+v)  aB  YY  aB  1-v  ^zz  oB  ~ 1-v  aB 

Upon  substitution  of  (26)  into  (5)  and  application  of  (4),  we 
readily  obtain 


M v6  h/2 

"as  ■ ^ 

-h/2 

A similar  calculation  reveals  that  ' 

_ N 6 - v6  - h/2 

"aB  ■ ^ 

Oa  - "’“"Oa 


(26) 


(27a) 


(27b) 

(27c) 


where 


in  “ m * l|/  “ 

aa  ^^0,0 


n » ii  « u 

oa  a, a 


h/2 

M ■ OE  J zTdZ 

-h/2 

h/2 

N„  » OE  / Tdz 

-h/2 


We  now  drop  (neglect)  the  integrals  in  (27a)  and  (27b) . Thus  we 
obtain  the  plate  stress-strain  relations: 


”o6  ' (i_v^)  " (1-v)  **o0 


- D(l-v)«^g  + vm6^3l  - 6^3 


Q„  - Ghq^ 


Upon  substitution  of  (6)  into  (29) , we  obtain  the  plate  stress- 
displacement  relations 


V ■ <"a.6  * “6.a>  * ^ "“*> 


"as  ■ '♦a.S*  *6, a>  * ^'^y./as'  ’ 


- « ffl><»a  * -,a> 


(30b) 


(30c) 
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Further  substitution  of  (30)  into  (21)  results  in  the  displacement 
equations  of  motion  for  the  plate: 


Eh  ‘"t’  O 


(31a) 


phw  - K Gh(l|-^^^  + W,^^)  + P 


(31b) 

(31c) 


Equation  (30)  can  also  be  written  in  (invariant)  vector  form: 


phu 


Eh 


2(l-v^) 


2^  ± -r  -*■ 

((l-v)V  u+(l+v)V(V*u)]  -- 


(1-v) 


(31a) 


ph^Vp  - [(l-vjV^tP  + (l+v)^(^*ili)]  - tc^Gh(il»+^w) 

(1-V) 


(31b) 


Phw  - K^Gh(^*lP  + 7^w)  + p 


(31c) 
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2.4  Energy  Considerations 

In  general,  the  time  rate  of  change  of  strain  energy  density  Is 
given  by  (see  p.  267  of  ref.  (11).  For  the  plate,  this 

can  be  expressed  as  p ( p ^ ‘ Consequently,  the 

time  rate  of  change  of  strain  energy  per  unit  area  of  plate  Is 


h/2 

■ .Jj  >W“aB*"ae’  * 


•«  M-mo  + Qq  + Nrtno 

oB  oB  oB  aB 


(32) 


«diere  we  have  used  (5).  However, 


9*r 

BiJ" 

B'if 

^p- 

“ Q •* 

aB 

E 

^‘*a 

3n^p  "aB 

combining  (32) 

and 

(33),  we  obtain 

fBtr 

V 

/ bU 

- »*aBJ 

>*aB 

M 

\ / 

- 2a)^ 

aB 

(33) 


This  result  must  be  valid  for  all  processes,  l.e.,  &nd  n^^ 

are  arbitrary.  We  thus  conclude  that 

3M  BU 


3*^ 
— E 


®a"  3a  ' 


N 


aB  3n 


(35) 


aB 


and  It  Is  evident  that  the  scalar  function 

- ftrp(mj^j,sij^2'*22'‘*l'‘*2''‘ll''*12'''22^  completely  characterises  the 
SMchanlcal  configuration  of  the  hosngeneous  plate.  In  view  of  (35) 
and  (29) , the  pertinent  strain  energy  density  function  for  the  plate 
is  given  by 
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2N  2 

2^/"  “ D(l-v)ni  -in  - + DVn  — — — m + n o n o + 

P aB  aB  (1-v)  (l+v)  aB  aB 


2n^  2 

71^  n + K Ghq^q^ 


(36) 


With  reference  to  (36)  and  (29),  we  also  note  that 
(>%/  * M 6in  o + Q 6q  + N -6n  _ 

p oB  aB  a ^ aB  aB 
and  this  result  was  used  in  (15). 


2. 5 Reduction  to  Classical  Plate  Theory. 

1 3" 

If  we  neglect  rotatory  inertia  in  (21),  i.e.,  if  — ph 
then  the  stress  equations  of  motion  reduce  to 


■ "aB.S 


or 


*'aB,Ba^  ‘ 


(38) 


In  classical  plate  theory  it  is  assumed  that  material  line  elements 
originally  normal  to  the  median  surface  remain  normal  during  deformation, 
i.e.,  Consequently  m^^-w^gand  i.e.,  the  plate  is  rigid 

with  respect  to  shearing  deformations.  In  this  case 


"aB 


aB 


(l-v) 


(39.) 


■ "(rt.B  ■ 


(39b) 


♦ ’ 

» or 

Q - -DV(V^w)  - (39b) 

Upon  substitution  of  (39)  into  (38),  we  obtain 

» 

phw  + dV^w  » p - — (40) 
Equation  (40)  is  the  well  known  equation  of  a plate  within  the  frame- 

> 

work  of  classical  plate  theory,  including  thermal  (bending)  terms 
(see  I2J,  p.  384,  eq.  (12.2,15)). 

I 

\ 

r 

f 

t 

r 

9 


t 
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3.0  RESOLUTION  OF  BOUNDARY  VALUE  PROBLEMS. 

3.1  Generalities 

Inspection  of  equations  (31)  euid  (22)  reveals  that  the  plate 
stretching  problem  can  be  decoupled  from  the  plate  bending  problem 
in  the  present  formulation.  In  the  present  section  we  shall  present 
a formal  solution  for  each  of  these  problems.  We  shall  use  the 
method  of  Williams  (a  modified  eigenfunction  technique)  to  obtain 
a solution.  Since  the  present  approach  is  a generalization  of  a 
technique  developed  in  reference  (1] , we  shall  present  only  the  final 
results,  omitting  most  of  the  details. 

3.2  Plane  Extensional  Motion. 

The  problem  is  characterized  by  (21c),  (22b),  and  (30a),  i.e.. 


where 


and 


Eh 


N_6 


2(1-V^) 


I(1-V)(U_  „+u-  ^)+2Vu..  6 


T aB 


'a,B  8,a''‘-''“Y,YaB'  (1-v) 


(u  or  N ) and  (u.  or  N .) 
n nn  x nx 


(41a) 


(41b) 


(41c) 


are  specified  on  C.  For  a unique  solution,  we  adjoin  the  Initial 
conditions 


The  solution  of  the  problem  characterized  in  (41)  is  given  by 

00 

a 1 2 o 12  a 12  1 


(41d) 


(42a) 


I 


i 

4 

I 
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J 


1*1 


.(i) 


(42b) 


The  terms  of  the  solution  (42)  which  carry  the  superscript  (s)  are 
obtained  from  the  quasi-static  problem 


«(s)  .0 
”a6,e  ° 


,(s) 


in  A 


N_6 


2(l-v^) 


^ [(1-v)(u‘^X®’)+2vu‘«>6  ^1  - 

a,B  B,a'  y,y  aB  (l-v) 


(43a) 


(43b) 


and 


(u^®^  or  N^®^)  and  (Uj*^  and  are 

n nn  x ^ x 


(43c) 


specified  on  C as  in  (41c) 

Quantities  in  (42)  which  carry  the  superscript  (i)  are  associated 
with  the  eigenfunctions  characterized  by  the  homogeneous  equations 


aB«B  i a 


in  A 


(u^^^  or  and  (ui^^  or  N^g^-O)  on  C 

n nn  x nx 


(44a) 

(44b) 

(44c) 


The  eigenfunctions  satisfy  the  orthonormality  relation 


[ “ '^ij 


With  the  above  definition,  it  can  be  shorn  that 


(44d) 


q^(t)  - (q^{o)-Q^(o)l  cosu^t  + jj-  (q^(o)-Q^(o)l  sinu^t 


+ / Qj^(T)  8inu)^(t-T)dT 


(45) 


I 


I 


I 
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where 


i i i (l-v)  o,a  ^nnnnnn£nS.£.n£. 

A c 

(46a) 


and 


q. {o)“Q.(o)  “ / phu^^^u^°^dA 
i 1 ' a ot 


qi(o)-Qi(o)  » / phu^^^u^°^dA 


(46b) 


(46c) 


3.3  Bending  (Flexural)  Motion. 

In  this  case  the  problem  is  stated  by  (21a),  (21b),  (22a), 
(30b)  and  (30c),  i.e.. 


-phw  + 0 + p = 0 

a, a 

1 3" 

■TTPh'l'  + Moo“Q  “0 
12  aB,6  a 


where 


and 


Oa  ■ " “'♦a  * “.o’ 

".6-  ?■>  ‘“-''><*a,8^6.a>«'^r,YV 


(ij»  or  M ),  (tj/j  or  M ) and  (w  or  Q ) 
n nn  x nx  n 


(47a) 

(47b) 

(47c) 

(47d) 

(47e) 


are  specified  on  C.  For  a unique  solution,  we  adjoin  the  initial 
conditions 


w(Xi,X2,0)  - w*°^  3 “ '^o°* 

w(x^,X2,0)  - w^°^  3 ^^(Xj^,X2,0)  - 


(47f) 


)d& 


The  solution  of  the  problem  characterized  by  (47)  is  given  by 
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(Xj^rX^.t) 

00 

- I 

i*l 

00 

- 1 

i=l 

VS 

8^”(Xl.rj.t) 

00 

+ 1 
i-1 

o 

00 

- 1 
i=l 

”a3 

(t) 


(48a) 


{48b) 


(48c) 


(48d) 


The  terms  of  the  solution  (48)  which  carry  the  superscript  (s)  are 
obtained  from  the  quasi  static  problem 


+ o - 0 
a ^ 

”a6,3  ° 


in  A 


where 


q‘®>  - .^Gh(,^‘®>tw<®>) 
*01  / ot 


and 


-rt’  - 5 (l-v)  “a6 


or  m'*’),  (♦!'>  or  m‘">)  „d  (»'•>  or 

n nn  x.  nt  n 


(49a) 

(49b) 

(49c) 

{49d) 

{49e) 


are  specified  on  C as  in  (47e) . Quantities  in  (48)  which  ceurry  the 
superscript  (i)  are  associated  with  the  eigenfunctions  characterized 
by  the  homogeneous  equations 


.(i) 


2.  (i) 


or-  - -pha,^w 


„(i)  .(1)  „ h"  ,.2  .(i) 

--Pl2“i*a 


’oB.e  '"a 


in  A 


(50a) 


(50b) 
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or  or  and  or  m‘^*=0)  (50c) 

n n nn  £ n£ 


on  C. 


The  eigenfunctions  satisfy  the  orthonormality  relations 


In  this  case  it  can  be  shown  that  q^(t)  is  again  defined  by  (45), 
and  in  this  case 


2n  ,4.x  L , <s)^(i)  (i),,(s)^,  (s)„(i)  ,(i)„(s) 

iJ.Q,  (t)  * 6 (w  Q -w  Q +ip  M -llf  M 
1 i ^ n n nn  ^n  nn 

- i ^ 7^  c>- 


q^(0)-Q.(0)  - / (phw<°’w‘^^  + P 12 


q^(o)-Q.(o)  - / + p 


In  the  following  we  shall  consider  the  special  case  of  a 
circular  plate  with  radius  a which  is  rigidly  clamed  at  its  boundary. 

The  plate  will  be  subjected  to  laser  irradiation  at  its  center,  i.e., 
the  excitation  as  trell  as  the  response  of  the  plate  will  be  axisymnetric. 
In  addition,  it  will  be  assumed  that  the  plate  is  at  rest  and  in  a state 
of  static  equilibrium  at  t>0.  For  these  reasons  we  now  vrrite  the 
principal  equation  in  axisymnetric  polar  coordinates  with  the  pole 
at  the  center  of  the  plate.  Consequently  we  have 
(a)  Plate  Stretching: 


u^  = uir  ,t) , U0  = 0 


(52a) 


! 


21 


Sn 

-Phu  + ^ - (N^^-Ngg)  = 0 


(52b) 


Eh  r 3u 


rr 


'ee 


l-v 


Eh 

l-v"* 


^ du  V » T 

^9F  ^ F - I^ 


3u  I . 

( V 37  + 7 - 


(52c) 


u(a,t)  = 0 

u(r,o)  = 0,  u(r,o)  = 0 


(52d) 

(52e) 


u(r,t)  » u*®^(r,t)  + I u'^'(r)*q  (t) 

i=l 

00 

N (r,t)  = Nj®Nr,t)  + I n^^Nr)-q.(t) 


(i) 


i=l 


Nee(r,t)  - NQ^gNr.t)  + I NgJ' (r)*q^(t)  ' 

i=l 


(i) 


2TT 


/ (ru'^^dr 


^ "T  3r 
o 

t 

qi(t)  - Q^(t)-(i»^  / Q^(T)sin  U)^(t-T)dT 


(53) 


(54) 


(b)  Plate  Bending 


9!^  8, 


-phw  + -r — + » 0 

3r  r 

- ^ ^ * r <'Vr-*%0’“2r 

-o(lit.^.)-^ 


•Vr  ■ “ ' ar  r 

•^e  - + F*  > ■ T 


0^  - K^Gh(V'  + 


T_ 

-V 


(55a) 


(55b) 
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4.0  LASER  INDUCED  HEATING  OF  THE  PLATE 

4.1  Ten^erature  Distribution  in  the  Plate 

In  order  to  obtain  a mathematical  model  of  the  temperature 

distribution  induced  in  the  plate  by  a pulsed  laser,  several 

assun^tions  will  be  made.  These  assumptions  yield  a model  which  is 

both  physically  realistic  and  mathematically  tractedile  in  the  forth- 

f coming  analysis.  The  geometries  of  the  clamped  circular  plate  and 

of  the  laser  irradiance  are  shown  in  Fig.  1.  With  reference  to  the 

dimensions  shown  in  Fig.  1,  the  five  basic  assumptions  are  as  follows: 

f (A-1)  The  laser  irradiance  is  specially  cucisymmetric  with  a Gaussian 

2 

distribution  of  the  form  exp[-(r/d)  ] where  d,  which  is  called 
the  laser  spot  size,  is  the  radial  distance  at  which  the 

* incident  irradiance  has  decreased  to  1/e  of  its  maximum  value  I . 

m 

(A-2)  The  laser  irradiance  also  has  a Gaussian  distribution  in  time 

2 

of  the  form  exp[-(t/T)  ) where  T,  which  is  called  the  pulse 
' duration,  is  the  time  required  for  the  incident  irradiance  to 

decrease  to  1/e  of  its  initial  value. 

(A-3)  The  radius  of  the  plate  is  essentially  infinite  compared  to 

* the  laser  spot  size,  i.e.,  a/d»l. 

(A-4)  The  plate  is  sufficiently  thick  so  that,  for  the  time  interval 

of  interest,  the  temperature  rise  of  the  back  surface  does  not 

^ have  an  appreciable  effect  on  the  notion.  This  condition  is 

difficult  to  express  analytically,  however  one  rough  approxi- 

2 

nation  is  obtained  fztm  the  Inequality  h /4Kt^>l.  Here  h is 

0 the  plate  thickness,  K is  the  thensal  diffusivity  of  the  plate 

and  t is  a characteristic  time,  e.g.,  the  fundaswntal  period 
P 


II 
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for  free  vibration  of  the  plate. 

(A-5)  Over  the  time  interval  of  Interest,  thermal  diffusion  in  the 

radial  direction  is  negligible.  This  condition  can  be  expressed 

2 

analytically  as  follows:  d /4Kt  >1. 

P 

Note  that  assun^tions  (A-3)  and  (A-4)  permit  us  to  treat  the  plate  as 
a semi-infinite  solid  when  computing  the  temperature  distribution. 

The  temperature  distribution  induced  in  a semi-infinite  solid  by  an 
incident,  Gaussian  laser  pulse  has  been  determined  by  J.H.  Bechtel  [3] . 
As  a result  of  assun^tion  (A-5),  Bechtel  shows  that  the  temperature 
distribution  is  approximated  by  the  function 


T t 

T(r,z,t)  = — exp[-(r/d)^]  / exp[-(tVT)^l 

rr 

2 ^ -1/2 
exp[-z  /4K(t-t') ] (t-t  ) dt 


(58) 


where  z is  the  distance  measured  into  the  solid  from  the  irradiated 
surface  and , 

T - I (1-R)  (TiK  (59) 

o m c 

In  (59),  I is  the  nuucimum  incident  Irradiance,  R is  the  reflectivity 
in 

of  the  plate  surface,  K , P,  and  C ^ure  the  plate's  thermal  conductivity, 

c 

density,  amd  specific  heat,  respectively.  Note  that  the  thermal 

diffusivity  K>K  /PC. 

c 

A sore  convenient  form  of  (58)  is  obtained  by  making  the  change 
of  variable  0-(t-t')/T.  The  result  is, 

OD 

T(r,z,t)  -T  exp(-(r/d)^)  / exp(-(7- 0)^)  •exp(-zV^KT6J0  ^^^d0 

O T 

O 


(60) 
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4.2  Thermoelastic  Forcing  Functions 

The  thermoelastic  forcing  functions  and  are  defined  in  (28) 
and  revrritten  here  as  follows: 


h/2  h 

N (r,t)  ■ OE  / TdZ  » OE  / Tdz 
^ -h/2  o 

h/2  h h 

M (r,t)  = aE  / TZdZ  - aE[/  Tzdz  - 7 / Tdz] 
-h/2  o ^ o 


(61a) 


(61b) 


where  Zoz-h/2  is  measured  from  the  middle  surface  of  the  plate  as 
shown  in  Fig.l.  Substitution  of  (60)  into  (61)  yields  the  following 
results: 

^ (62a) 

(62b) 

X u 

where 


N (r,t)  * N exp[-(r/d)  )N(t) 
T O 

M^(r,t)  exp(- (r/d)^)M(t) 


1/2 

N « TT  ' OEI  (1-R)T/PC 
o m 

M » N h/2 
o o 


(63) 


N(t)  / exp [-(7- 6)^1  erf  (<|))de 

o 

00 

M(t)  / exp[-(7-e)^){erf(())  )-2n'^''^(l-exp(-<|)^)l/«j)}de 


(64a) 


4»  - h(4KTe) 


o 
-1/2 


(64b) 

(64c) 


4.3  Numerical  Values  For  Parameters 

In  order  to  coiq>are  the  forthcoming  analytical  results  for  the 
motion  of  the  plate  with  the  experimental  results,  the  following 
numerical  values  will  be  used  for  the  parameters: 


O 
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j E - 46.7  X 10^  Ib/in^ 

V » 0.33 
2 

K » 0.86 

p - .7287606  X lo“^  Ib'sec^in^ 

K « .726  X 10  ^ in^/sec. 

h - .0195  in. 

f a > 4.523  in. 

A “ a/d  » 16 
-9 

T • 40  X 10  sec. 

^ These  parameters  correspond  to  the  experimental  phase  of  this  project, 

as  described  in  Section  7 of  this  report. 

I 


I 

i 

1 

F 

f ^ 


28 


5.0  FREE  VIBRATION  OF  CLAMPED  CIRCULAR  PIATES. 


5.1  Transverse  Vibrations 


In  the  case  of  axisymmetric , free,  transverse  vibrations, 


where  R'^r/a  is  a dimensionless  radial  coordinate.  The  stress- 
displacement  relations  become  (see  (55b)): 


Substitution  of  (66)  into  (55a)  yields  the  following  homogeneous 
system  of  equations  for  the  determination  of  the  eigenfunctions: 
For  0SR<1, 


(67a) 


R dR 


(Oj  - 

^ R 


a dR 


while  at  R>1, 


(1)  - - 0 (67c) 

In  addition  «ra  require  w^^^  (0)  and  i|'^^No)  to  remain  ]x)unded.  In  (67), 


tc^G(l-v^)  K^(l-V) 


(eeb) 


fi 


(68c) 


2 

V 

P 


E 

P(l-V  ) 


The  normalization  condition  (50)  Isecomes: 


(68d) 


2Tia^ph  / ((w^^^)^+H^a^(V^^)^)RdR  - 1 (69) 

o 

The  solution  of  (67)  through  (69)  is  given  in  [1]  (pp.  483-488). 

The  following  is  a sumnary  of  that  solution  in  a form  amenable  to 
numerical  computations. 

When  solving  equations  (67),  one  naturally  obtains  a division 

of  the  frequency  spectrum  into  two  intervals,  the  low  frequency 

2 2 2 2 
range  (0<f.H  Q <1) , and  the  high  frequency  range  (J.H  ft  >1) . The 

number  of  natural  frequencies  contained  in  the  low  frequency  range 

increases  as  the  thicliness  of  the  plate  decreases.  In  fact,  in  the 

forthcoedng  forced  motion  example  we  will  find  that  all  the  natural 

frequencies  which  contribute  significantly  to  the  solution  are 

contained  in  the  low  frequency  range.  Therefore,  the  forthcoming 

analysis  will  be  restricted  to  the  low  frequency  range 

0 < iH^n^  < 1 (70) 


Tlie  frequency  equation  as  well  as  the  eigenfunctions  can  be  expressed 
conveniently  in  terms  of  the  parameter  X «rhich  is  related  to  the 
frequency  ft  as  follows:  Let, 


C - 1 ♦ (t+l)X 


(71a) 

(71b) 


r 
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r 

r M * 2£aVc^  (71c) 

Then, 

j ) (71d) 

♦ 

In  the  low  frequency  range,  y«l  and  therefore  it  is  advisable  to  use 
a truncated  Taylor  series  to  compute  the  values  of  /l-2u.  Substitution 
of  this  series  into  (71d)  yields  the  alternate  formula 

r 

- A^n  (72a) 

where , 

r n •=  (72b) 


and, 

C ~ l+y(.5+y (.5+p(.625+W(.875  + U(1.3125  + y (2. 0625) ) ) ) ) ) (72c) 

In  (72c),  a sufficient  number  of  terms  have  been  retained  in  the 


Taylor  series  to  yield  values  of  ^ correct  to  at  least  eight 
significant  figures.  In  view  of  (72a)  and  (71a)  we  conclude  that 

fi  - A (73) 

H 

This  last  equation  is  very  revealing.  For  the  lowest  frequencies  of 
a thin  plate,  A«l,  therefore  and  O-HX  . 

In  terms  of  the  parameter  X,  the  frequency  equation  may  be 
written  as  follows i 

J (X)I, (Y)  + f(X,Y)J-{X)I  (Y)  - 0 (74) 

o I 1 o 


In  (74),  Y Is  related  to  X as  follows: 


bet. 


c 


1-fl 


(75a) 
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P - 1 + (75b) 

and, 

6 = -2e(l-e)/p^  (75c) 

Then, 

2 

e(|)  - jp  (/l-26  -1)  (75d) 

Since  6«1  in  the  low  frequency  range,  /l-26  may  be  computed  from  a 
truncated  Taylor  series,  in  which  case  (75d)  is  replaced  by  the 
following  formula: 
y 2 

(— ) = r = (l-e)q/p  (76a) 

where , 

g S 1+6(.5+6(.5+6(.625+6(. 875+6(1, 3125+6(2. 0625))))))  (76b) 

Consequently, 

y«sX  where  s = (77) 

Note  that  for  the  lowest  frequencies  of  a thin  plate,  s=l  and  therefore 
y>X.  The  function  f(X,y)  appearing  in  (74)  is  defined  as  follows: 

f(X,y)  - s/b  (78a) 

where 

b - (r+A)/(l+r«,)  (78b) 

Note  that  f(X,y)>l  for  the  lowest  frequencies  of  a thin  plate.  Com- 
bining this  with  the  observation  that  y>X,  the  lowest  frequencies  of 
a thin  plate  nay  be  approximated  by  the  roots  of  the  equation 

J (X)I- (X)  + J, (X)I  (X)  - 0 
o 1 1 o 

which  is  in  fact  the  frequency  equation  of  classical  plate  theory  for 
a claaq;)ed  circular  plate. 

y 

He  note  that  the  functions  e , X (Y) , and  I (Y)  all  Increase 

O X 
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exponentially  for  large  values  of  Y and  have  no  positive  real  roots. 

Therefore,  the  frequency  equation  (74)  can  be  modified  by  dividing 

-y 

the  entire  equation  by  e I (Y)!  (Y)  without  affecting  its  roots. 

o 1 

This  yields  the  modified  frequency  equation 
J (X)  J (X) 

F(X)  - — + f (X,Y)  - 0 (79) 

e I (Y)  e"  I,  (Y) 

o 1 

Equation  (79)  is  more  convenient  to  use  in  numerical  computation  than 
(74)  since  the  maximum  amplitudes  of  the  oscillations  in  F(X)  remain 
bounded  as  X increases.  Furthermore,  an  examination  of  the  asymptotic 
behavior  of  Bessel's  functions  reveals  that  the  roots  of  (79)  are 
approximately  given  by 


X^  » iTT  for  i»l,  2,  3,  ... 

X.  S (i-'7)Tr  for  i=N,  N-1,  N-2,  ... 

X 4 


(80) 


where  N represents  the  total  number  of  frequencies  in  the  low  frequency 

2 2 

range.  The  value  of  N is  easily  computed  by  noting  that  £h  (I  <1  implies 

that  e<l  which  in  turn  Implies  that  X < ^ /(£,+l)/S,.  Therefore,  N is 

the  largest  integer  which  is  less  than  ^ »^(i+l)/A.  Using  the 

parameters  given  in  4.3  yields  N>290  in  the  present  example.  The 

first  35  roots  of  (79)  were  ccm^ted  using  a modified  Ne%/ton-Raphson 

Iteration  technique  with  starting  values  of  X^,...x^g  given  by  the  first 

equation  of  (80).  The  modification  «ias  to  use  the  first  central 

difference  [F (x-f. 0001) -F(x-. 0001 ]/. 0002  to  approximate  the  derivative 
dP 

^ in  the  iteration  formula.  This  technique  yielded  the  roots  to  at 

least  eight  significant  figures  accuracy  within  t%«o  or  three  iterations. 

“i 

The  first  35  values  of  X. , , and  — are  listed  in  Table  I. 
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The  corresponding  normalized  eigenfunctions  are  given  by  the  following: 
For  OiRil, 


(81a) 

- X.  — (r) 

1 a 1 

(81b) 

- (-^)hH\.  (R)-(1-V)X.F.^^Nr)/R1 

" iV  ^ 1 

(81c) 

» {-^)hH\^I\«jF^^^R)  + (l-V)X^F|^^R)/R] 

(81d) 

1-V  2 

— (-^)  -X.  jpFf^^R) 
r , .2  a 1 *.  X 

1-V  1 

(81e) 

where , 

= J (X  R)-c.I  (Y.R) 

1 o i 1 o X 

(82a) 

t 

F?^^  * a. IJ  (X  R)-d  I, (y.R)l 

1 ill  i 1 i 

(82b) 

^ '1<V' 

(82c) 

r 

'1"’  ■ 

(82d) 

and. 

2 1 

21..=pwf 

(83a) 

f 

(83b) 

f 

1 2 »l 

*i  • ‘ ? \*1  " 

(83c) 

■i  - * v! 

(83d) 

n 

2(l-r  f ) 2 2 

<^1  ■ «^(l«;)  ll*(l«)*i\n,X^/,2r^)) 

(83e) 
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a.  - (1+r.i,)/ (!+«,) 
1 1 

b.  - (r.+Z)/(l+r  £) 

1 1 i 

c.  - J (X.)/I  (YJ 
1 o 1 o i 

d.  = J,  (X.)/I,  (YJ 

1 1 1 1 i 


6 = s d. 

1 11 


o ,7r  2 

“i  ° ^2 


-1 


fi  = f(X..Y.) 


5. 2 Extensional  Vibrations 

In  the  case  of  axisymmetric,  free,  extensional  vibrations. 


<i>  - « 

U0  = 0 


rr 


,(i) 

00 


Eh 


a(l-v") 

Eh 

a(l-v^) 


2.  dR  R ' 


' dR  R ’ 


Substitution  of  (86)  into  (52b)  yields  the  following  homogeneous 
system  of  equations  for  the  determination  of  the  eigenfunctions: 
For  01R<1, 


while  at  R-1, 

..(i) 


(1)  - 0 


(84) 


(85) 


(86a) 


(86b) 


(87a) 


(87b) 
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We  also  require  u'  ' {0)  to  remain  bounded.  The  corresponding  normali- 
zation condition  (44d)  becomes, 


2TTa^ph  / (u^^^)^RdR  = 1 


The  solution  of  (87)  ^md  (88)  can  easily  be  shown  to  be 


u (R)  = A^Jj(n^R) 


(R)  = ^ [fi.J  (fi.R)  - J,  (fi.R)] 

...  EhA.  . 

ad-v^)  ^ ° ^ ^ ^ 


where , 


“ Tia^phJ^Cn.) 
2 o 1 


J^(n^)  - 0, 


i^ 1,2,3,...  . 


An  examination  of  the  asymptotic  form  of  Jj^(R)  reveals  that  the 
natural  frequencies  are  approximately  given  by 


n,  « (i  + ' 

i 4 


i-1,2,3,, 


Table  II  lists  the  first  30  natural  frequencies  ta)cen  from  [4]  (p.  748 

Table  VII).  The  corresponding  values  of  J (0  ) and  (£).)! 

ol  iziol 

are  also  listed  for  convenience  in  future  computations. 
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f 

I 


TABLE  II 

NATURAL  FREQUENCIES  OF  EXTENSIONAL  VIBRATION 


1 

J (fi.) 

O 1 

a? 

1 

1 

3.8317060 

-.4027593957 

1.024227855 

2 

7.0155867 

.3001157525 

1.007484896 

3 

10.1734681 

-.2497048771 

1.003591665 

4 

13.3236919 

.2183594072 

1.002101537 

5 

16.4706301 

-.1964653715 

1.001377622 

6 

19.6158585 

.1800633753 

1.000972230 

7 

22.7600844 

-.1671846005 

1.000722608 

8 

25.9036721 

.1567249863 

1.000558091 

9 

29.0468285 

-.1480111100 

1.000443971 

10 

32.1896799 

.1406057982 

1.000361582 

11 

35.3323076 

-.1342112403 

1.000300165 

12 

38.4747662 

.1286166221 

1.000253167 

13 

41.6170942 

-.1236679608 

1.000216398 

14 

44.7593190 

.1192498120 

1.000187095 

15 

47.9014609 

-.1152736941 

1.000163364 

16 

51.0435352 

.1116704969 

1.000143877 

17 

54.1855536 

-.1083853489 

1.000127682 

18 

57.3275254 

.1053740554 

1.000114073 

19 

60.4694578 

-.1026005671 

1.000102530 

20 

63.6113567 

.1000351468 

1.000092653 

21 

66.7532267 

-.0976530158 

1.000084139 

22 

69.8950718 

.0954333390 

1.000076746 

23 

73.0368952 

-.0933584533 

1.000070287 

24 

76.1786996 

.0914132722 

1.000064609 

25 

79.3204872 

-.0895848220 

1.000059593 

26 

82.4622599 

.0878618760 

1.000055139 

27 

85.6040194 

-.0862346634 

1.000051167 

28 

88.7457671 

.0846946348 

1.000047609 

29 

91.8875043 

-.0832342730 

1.000044408 

30 

95.0292318 

.0818469379 



1.000041521 

i 
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6.0  FORCED  MOTION  OF  THE  CLAMPED  CIRCULAR  PLATE 


6.1  Transverse  Motion 


In  the  subsequent  forced  motion  analysis  we  shall  find  it  con- 
venient to  introduce  a dimensionless  time  scale  based  on  the  fundamental 
period  of  free  transverse  vibrations.  Thus,  we  define 


t,  * t;:  t - (119.96775)t 
1 2tt 


Ola) 


or,  conversely. 


— t “ (8.3355736  x 10"^)t, 
1 1 


(91b) 


The  thermoelastic  bending  moment  M^  is  given  by  (62b)  with  M(t) 


defined  in  (64).  A very  accurate  closed  form  approximation  of  M(t) 


t 2 

ciui  be  obtained  from  (64b)  as  follows.  The  term  expl-(— -6)  ) tends 


to  strongly  focus  our  attention  on  those  values  of  6 for  which 

.-11 


|^-0l<5.  Indeed,  when  |^-e|i5  then  expI-(^- e)^Jil.4  x 10 
and  the  contrDsutions  to  the  Integral  in  this  range  are  negligible. 


The  smallest  time  of  interest  in  the  forced  motion  solution  is  on  the 


order  of  of  the  fundamental  period.  Thus,  the  time  interval  of 


Interest  corresponds  to  t.i.Ol  or  t28.336  x 10~^  or  - i 2084.  At  —=  2084, 

1 T T 


the  lnequalj.ty  |^-0|<5  yields  the  inequality  2O79<0<2O89  which  when 
substituted  into  (64c)  yields  the  result  12.52<(^<12.54.  Note  that  ^ 
is  essentially  constant  over  the  range  of  values  of  6 required  to 


evaluate  the  integral^ (64b)  for  this  instant  of  time.  Furthermore,  if 


— > 2064  then  ^ varies  even  less  over  the  range  of  values  of  0 needed  to 


perform  the  integration.  Consequently,  for  t^^.Ol,  ^ can  l>e  considered 


a constant  in  the  Integral  (64b)  with  9 • i.e. 
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-1/2 


or 


4>  » 1.2533390  t. 


-1/2 


(92a) 


(92b) 


After  removing  the  constant  expression  involving  ((i  from  the  integral, 
all  that  remains  to  be  integrated  is 

OD 

-1/2 


It  ^ / expI-(^-S)  Jd0  * ■^(l+erf(^)J  * 1 


since  — i 2084  and  erf(x)  is  one  to  ten  decimal  places  for  any  x>5. 
We  thus  arrive  at  the  following  approximation  for  the  thennoelastic 
bending  moment: 


M{tj)  - erf  {<|»-2Tr"^''^[l-exp(-4)^)]/4> 


(93) 


According  to  (64b),  M(0)*l/2  but  M(t)  very  quickly  increases  to  1. 
Since  this  rise  time  is  less  th£m  of  the  fundamental  period  we 
will  simply  assume  that  M(0)>'1  as  predicted  by  (93).  As  a result  of 
these  observations  we  shall  henceforth  use  (62b)  with  M(tj^)  given 
by  (93)  to  represent  the  thermoelastic  bending  moment.  A graph  of 
M(tj^)  vs.  tj^  is  shown  in  Pig.  (2). 

The  generalized  force  Qj^(tj^)  is  now  evaluated  by  substituting 
(62b)  into  (57).  The  result  is, 

o 


. H(t . / 

1“V  O 1 ' 


(94) 


where  A>^  is  the  spot  ratio.  Consider  the  integral 


<►  % 


O 
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I . / 


(95) 


By  substituting  (81b)  and  (82b)  into  (95)  we  obtain 
1 2 2 

I - ^ Vi“i  ^ ^ IJ^(X^R)-d^I^(y^R)JdR 

O 

2 2 

R*^  — •»!  1 O 

Notice  that  at  R“l,  e *e  *6.6  x 10  and  therefore  the  upper 
limit  of  the  integral  can  t>e  extended  to  infinity  without  noticably 
affecting  the  value  of  I.  This  technique  allows  us  to  obtain  the 
following  closed  form  expression  for  I: 

I - ^ i - {X,  expt-(X  /2A)^)  - d.Y,  exp[(y,/2A)^])  (96) 

.,4  1 1 11  i 

4A  a 

Substitution  of  (96)  into  (94)  then  yields  the  desired  result, 


TTM  (1+V)  K 

hw  ■ "“i’ 

4A  Eh  H 


where 


- e^{exp(-(X^/2A)^J-6^  expI(Y^/2A)^]} 


(97a) 


(97b) 


with  e^  and  3^  defined  in  (84).  Further  substitution  of  (97)  into 
the  second  of  (57)  yields  the  result 
1TM  (l+v)  K, 


9j^(tj^) 

4A  Eh  H 


(98a) 


where, 


gi(ti)  ■ M(tj^)-2irR^  M(8)sin  2TrR^ (tj^-e)de 


(98b) 


and  where 


(4. 


(98c) 


An  efficient  schesM  for  the  numerical  evaluation  of  9^  (^2.^  given 


t 


in  Appendix  I. 

The  quasi-static  problem  for  the  present  .example  is  characterized 
by  the  following  equations: 


For  0 i R < 1, 

-0 

R dR  r 


dR  * R '”rr  ”00  ^ *°r 


where , 


.(s)  1 2_.  . . (s)  dw'“'  , 

«t  ■ r " * “s~  ’ 

(s) 

rr  a dR  R ^ 1-v 


„(s)  .D  1 ^(s), 

00  a dR  R ^ ' 1-v 


(99a) 


(99b) 


(100a) 


(100b) 


(100c) 


The  corresponding  boundary  conditions  are: 


w‘®^  (1)  - (1)  - 0 

while,  w^^^  (0)  a!!»3  i^^^No)  remain  bounded. 

The  solution  of  (99)  through  (101)  is  sunnarized  Ijelow. 


w^®^  (R.tj^) 


(R,t^) 


M a 

M(t,)W  (R) 

DA^l-V)  ^ « 

M a 

—=-2 M(t,)'l'  (R) 

DA^l-V)  ^ ® 


(102a) 


(102b) 


(102c) 


•*06  <*'V  “T”‘V”0S<"’ 


(102d) 


H 0 


(102e) 
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I » 

1 

r 


9 


A 


r 


f 


(I 


where , 


W (R)  - - 7[(1-R^) (l-e‘^  ) + InR^  + E, (A^R^)-E, (A^) J 
s 4 11 

2 2 2 

(R)  - )-R{l-e~^  )] 

S A K 

-h  " -h 

A R A 

2e"''  ^ d-e"^  ^ ) 

A R A 


(103a) 


(103b) 


(103c) 


(103d) 


In  (103a),  E^(X)  is  the  exponential  integral  defined  as  follows: 


oo  -t 

e 


Ei(X)  - / ^ dt  » / d-e"’')  ^ - Inx-Y 


(104) 


where  Y*. 5772156649  ...  is  Euler's  constant.  In  the  limit  as  R^O, 
(103a)  becomes 

12  2 -A^ 

N (0)  - ^IlnA  + E,  (A  ) + Y - d-e  ))  (105) 

S 4 i. 


Fig.  (3)  shows  the  radial  variation  of  the  q;uasi-static  deflection 

H (R) . As  expected  the  maximum  deflection  occurs  at  the  center  of  the 
s 

plate.  Pig.  (4)  shows  the  radial  variation  of  the  quasi-static  bending 
moments.  Again  the  siaximum  stoment  occurs  at  the  center  of  the  plate 
where  M . The  coBq>lete  forced  motion  solution  is  now  obtained  by 

RS  vS 

substituting  (102),  (96),  and  (81)  into  (48).  The  results  are 
M (1+v) 

w(R,t-)  - - lW_(R)M(t,)-S, (R,t,)l  (106a) 

^ A^H^Eh  ® 111 

N 

Mrr(R»t^)  - -f-  + S2(R,t^)l  (106b) 

M 

MeQ(R,ti)  - -y  lMgg(R)M(t^)  + Sj(R,t^)I  d06c) 

M 

“ • 2d(l=^  S^(R.t^)  (106d) 


I 


V = 0.33 


Figure  4:  Quasi  - Static  Bending  Moments  and  Normal  Forces 
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where , 


w 


(107a) 


TT  r 2 

=2<''-h’  ■ 7T  Vi'-ii^ O-  ”i'‘i> 

2A  l«l  1 

'^'"■'1’  ■ J-,  Vi'-i^r—  * ^?jr-'’i'S> 

2A  i»l  1 


(107b) 


(107c) 


oo  n 


i"'''i>  -rr  i,  ’‘i3r''”’"'”i''i> 


(107d) 


2A  i=l  H X. 

1 


Since  the  maximum  deflection  emd  moment  occurs  at  the  center  of  the 
plate,  we  set  R*0  in  (107)  to  obtain: 


S (O.t  ) = I (t  ) 

i-1 


(108a) 


SjfO.t  ) - I s‘^^g  (t  ) 
i-1 


(108b) 


Sj(0,ti).  - I S^^’g^(t^)  - 8^(0, t) 


where , 


S<i>  . J K^i) 


s<i)  - 2Lk  ^ (-2—- 

“2  2 1 ^2  ' 1-  V °3  ' 

„(i)  W , *i  ,2  ^ „(1),  „(1) 

®3  “ 2 *^i  ]T  ‘T^  * ^2  ^ " ®2 


(108c) 


(108d) 


(108e) 


(108f) 


and 
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1-c  . 


1 


(108g) 


- a^d+b^c^)  (108h) 

j a^(l-B^)  (108i) 


The  first  35  values  of  presented  in  Table  III. 

This  table  shows  that  by  suiraning  the  first  35  values  we  are  guaranteed 
at  least  4 significant  figures  of  acctiracy  in  our  final  results. 

Fig.  (5)  is  a graph  of  the  dimensionless,  negative  displacement  at  the 
center  of  the  plate: 


W(0,tj^)  = Wg(0)M(t^)  - Sj^(0,tj^) 


(109) 


The  negative  displacement  was  used  to  facilitate  later  cosQ>arison  with 
the  experimental  results.  The  periodic  nature  of  S^(0,t^)  is  evident 
in  this  graph  with  the  second  cycle  beginning  at  approximately 
t^>1.04.  The  large  initial  slope  of  the  deflection  curve  indicates 
that  the  thermoelastic  forcing  of  the  plate  behaves  somewhat  like  an 
inqinilse  (finite  initial  velocity-slope)  or  perhaps  even  a doublet 
(infinite  initial  velocity-slope).  Fig.  6 shows  the  dimensionless 
bending  moment  at  the  center  of  the  plate 


M 

rr 


(0,t^) 

M 


o 


M 


|(MRs(0)M(ti)  + S2(0,t^)] 


(110) 


The  very  nervous  behavior  of  the  bending  moment  reinforces  the  previous 
observation  concerning  the  singular  nature  of  the  thermoelastic  forcing 
of  the  plate. 
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TABLE  III 


MODAL  CONSTANTS  FOR  TRANSVERSE  MOTION 


1 

5.71  1228377“! 

4 .7355160’^“! 

3 . 04  56808F“l 

2.3860140F  1 

1 . 9 3 8 4 6 6 5 77  1 

1 . 5220  37  5F“l 

1 . 3612248P“l 

1.0690083'^  1 

9.9915163P  2 

7 . 8473202F“2 

7 . 401331077“2 

5.  8836776F“2 

5. 666174777"2 

4.it502032F  2 

4 . 2911681P  2 

3. 3702756F“2 

3.23R2838P  2 

2.5433422r  2 

2 . 427010277“2 

1 . 9061693F“2 

1.802270877  2 

1.4155002F  2 

1 . 3237398£’“2 

1 . 0396628F“2 

9.60381  6577“3 

7 . 54  2 8198F“3 

6. 875373477“3 

5 . 3999057F“3 

4 . 8529442/7“3 

3 . 8114g35F“3 

3 . 375087377“3 

2. 6507874F“3 

2. 311  541  677“3 

1. 8154805 F“ 3 

1 . 558338677“3 

1.2239162F“3 

1 . 033720977“3 

8.1188251F“4 

6.745111877  4 

5. 2975984F“4 

4 . 3281518F“4 

3 . 3993225F“4 

2.7305041F“4 

2 . 1445329F“4 

1.6932583F“4 

1 . 3298820F“4 

1.0319717F“4 

8.1050872F“5 

6. 1803085F“5 

4.8540029F  5 

3 . 6365677F“5 

2. 8561536F“5 

2. 1021345F“5 

1 . 6510126^“5 

1.1936320F“5 

9. 3747641F“6 

6. 6570270F“6 

5. 228416RF“6 

3 . 6463055F“6 

2.8630017F“6 

1 . 9613512F"6 

1 . 5404416F“6 

1 . 0359937F“6 

8.1366755F“7 

5.373175BF“7 

4 . 2200824F“7 

2.7362285F“7 

2. 1490289F“7 

1 . 3680389F“7 

1 .074455?F“7 

T .35505  2 577^2 
7.3037055  7:  2 
1.0512335P"! 
1. 3101330F"! 
1.5008  00  577  1 

1 . R1801Uii-“i 

1.0051482^  1 

1 . 04  5555: 
1.570072ir  1 

1 .4  5117  5277“! 
1 . 30240047’"! 
1 . 1369863.F“1 
9. 067807  0.>7  2 
8 . 01  531  9577"2 
0.4  84  59157’  2 


, 1 A z r ^ 

3.9538539F“2 
2. 982826477“2 
2. 20034077’“2 
1 . 5B70155F“2 
1 . 1 20760477“2 
7 . 7427371P“3 
5. 23  5754877“3 
3.4661286F"3 
2.246771077  3 
1 . 4262O520’"3 
8. 86681  9277“4 
5. 399632477"4 
3 . 221172477"4 
1.882590477  4 
1 .078015077“4 
6.04856R077“5 
3. 3255850P“5 
1 . 7918336P"5 
9.4616222P“6 
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A = 16 

t = (8.33557  x10‘’)t 
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6. 2 Extensional  Motion 

As  in  6.1,  we  shall  introduce  a dimensionless  time  scale  based  on 
the  fund2unental  period  of  free  extensional  vibrations,  i.e. , 

“l 

t_  = — t = (36156. 676)t  (111a) 

2 2tT 

or  conversely, 

t = ^ t^  = (.027657410  x 10“^)t2  (111b) 

An  accurate  approximation  of  the  thermoelastic  normal  force  can  be 

obtained  in  a manner  similar  to  that  used  in  section  6.1.  The  reasoning 

goes  as  follows.  The  time  interval  of  interest  in  the  present  problem 

corresponds  to  the  first  few  cycles  of  the  extensional  motion.  Thus, 

we  shall  consider  the  interval  .Olit  ilO  which  implies  6.915  i i 6915. 

2 T 

Since  l0-^|<5,  1.915^0^6920  and  therefore  6. 876i<|)i413 . In  this  interval, 
erf((^)'=l  to  more  than  ten  decimal  places  and  therefore  (64a)  becomes, 

OO 

N(t)  - j exp[-(Y-  0)^]d0  - |(l+erf(|))  = 1 

o 

since  erf  (^)  is  also  equal  to  one  in  the  time  interval  of  interest. 

Thus,  we  shall  use  N(t)»l  for  t>0,  or 


N^(R,t)  “ expl-(r/d)  ] for  t>0 


(112) 


The  generalized  force  obtained  by  substituting  (112) 

and  (89a)  into  (54) . The  result  is 

2 i * f ^ “A  R 


“i®i‘^2^  - RdR. 


2 2 

“An  -112 

Again,  noting  that  e «6.6  x 10  at  R*l,  we  may  extend  the  upper 


f 
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limit  in  the  above  integral  to  infinity  without  noticably  affecting  the 
value  of  the  Integral  to  obtain: 


/ J (ft.R)e"^  ^ HdR  i e ^ 
o ° " 2a2 


2 2 

-n:/4A^ 


Thus, 


N aitA  d fi.  2 

i P 


(113a) 


N anA  d n.  2 

‘ii<’^2>  ■ ^TTR.t^ 

i P 


{113b) 


vrtidTC f R , “ 9 i*l • 

i U)^ 

The  corresponding  quasi-static  problem  is  characterized  by  the 


equation 


J-rlA.  ,o„<s),,  ,1+V,  ^ 

dR  R dR  ^ ^ Eh  * dR  ' 


0:SR<1 


(114) 


(s ) (s ) 

where  u (1)>0  and  u (O)  is  Isounded.  The  solution  of  (114)  is  readily 
found  to  be 


(s)  N (1+v) 

u (R) 

d 2Eh  s 


(115a) 


tdiere. 


2 2 2 

U (R)  - AR(  d-e"^  d-e"^  )1 

® A^R^  A^ 


(115b) 


Substitution  of  (115)  into  (52c)  yields  the  normal  forces 


o 


(116a) 
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I cneb) 


where  Npg(R)*Mj^{R)  and  NQg(R)-MQg (R) , (see  eq.  (103)).  The  radial 
variation  of  the  dimensionless  displacement  U (R)  is  shovm  in  Fig.  (7). 
Note  that  the  maximum  displacement  occurs  approximately  at  R=0.07  and 
thus  that  is  where  we  shall  study  the  motion  of  the  plate. 

The  total  solution  is  now  obtained  by  substituting  (115) , (116) , 
(113),  2md  (89)  into  (53).  The  result  is: 


u(R,t^)  N (1+V) 
2 O 


[U^(R)-S, (R,t^)] 


2Eh  ‘“s'"'  “1'  '2 


1 S (R,t  ) S (R.t  ) 


N (R,t  ) , S (R,t  ) VS  (R,t  ) 


where , 


(«y2A)‘ 


Sj^(R,t2)  * ^ I O^Jj^(n^R)e  ^ cos  2ffR^t 


i-1 


“ -(fi  /2A) 

A i“l 


and 


(117a) 

(117b) 

(117c) 

(117d) 

(117e) 


«?  “ I?  fl,J^(n,)l“^  (117f) 

1 2 10  1 

The  dimensionless  displacement 

u(.07,  t^)  - 0^(.07)-Sj^(.07,  t^) 

is  shown  in  Fig.  (8) . The  wave  nature  of  the  solution  is  clearly  evident 
in  this  graph.  The  waves  originating  at  the  heated  spot  propagate  out  to 


= 0.07 


iKR.tj)  vs.  t 
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the  clamped  boundary  are  reflected  and  then  propagate  back  again 
towards  the  center  of  the  plate.  The  passage  of  these  ««aves  at 
R“.07  is  clearly  seen  in  this  graph.  Since  the  maximum  normal 


forces  occur  at  the  center  of  the  plate,  we  plot 


N (0,t_) 
rr  2 


“ee'o-'j' 


(118) 


in  Fig.  (9).  In  view  of  (117), 


(119a) 


(119b) 


(119c) 


(119d) 


An  examination  of  reveals  that  30  terms  are  sufficient  to  obtain 

S^(0,t^}  to  three  significant  figures  accuracy.  The  dimensionless 

normal  force  N (0,t^)  is  shown  in  Fig.  (9). 
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7,0  EXPERIMENTS. 

The  measurements  performed  during  the  initial  phase  of  the  study 
of  Laser  interaction  with  finite  structures  have  been  based  on  the 
concept  that  the  object  of  the  investigation  is  to  provide  a careful 
conparison  between  experimental  data  concerning  Laser  generated 
thermoelastic  defonritions  and  stresses  in  finite  structures  with  the 
theory  of  such  processes. 

It  is  clear  that  a considerable  quantity  of  variable  quality 
experimental  evidence  has  indicated  that  the  mechanical  impulse 
transferred  to  a target  by  a Laser  pulse  could,  in  general,  be 
understood  in  terms  of  material  evaporation  and  removal,  and  surface 
plasma  generation  with  and  without  subsequent  plasma  heating, 
possibly  involving  the  creation  of  LSC  g LSD  waves  in  the  plasma. 

Most  of  these  studies  have  involved  the  use  of  either  piezoelectric 
detectors  or  ballistic  impulse  detection  of  gross  target  motion, 
both  in  vacuum  and  in  air.  However,  at  least  in  the  unclassified 
literature,  relatively  little  effort  has  been  directed  toward  an 
understanding  of  the  Laser  excitation  of  elastic  oscillations  in 
finite  structures.  Since  these  phenomena  are  of  significance  in  the 
theoretical  treatment  of  plate  dynamics,  and  since  induced  stress 
waves  and  target  oscillations  could  be  of  practical  value,  we  felt 
that  a careful  comparison  of  theory  and  experiment  would  be  both 
useful  and  interesting. 

It  is  immediately  clear  that  any  study  of  the  dynamic  motion 
of  a finite  structure  could  involve  the  thermoelastic  induced  motion 
of  the  target  in  addition  to  any  motion  induced  via  a plasma  shock  or 
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material  blow  off.  Therefore,  our  approach  to  this  problem  has 

involved  a detailed  theoretical  and  experimental  Investigation  of 

the  purely  thermoelastic  induced  motion  of  a circular  thin  plate 

clanqped  in  such  a manner  that  the  plate  has  zero  displacement  and 

slope  at  its  outer  free  edge.  The  dimensions  of  the  clamp  ring  are 

outer  radius  14.6cm,  inner  radius  11.45cm,  thickness  5.7cm,  and 

mass  11.5kg.  This  clamp  ring  is  mounted  in  a horizontal  plane  by 

means  of  three  support  rods  which  serve  to  connect  the  clan^  mount 

to  a. horizontal  plate  loaded  uniformly  with  460kg  of  lead.  This 

technique  is  necessary  because  access  to  the  rear  surface  of  the 

target  plate  is  required,  and  the  Laser  inpulse  is  easily  cap2d3le 

0 

of  inducing  large  amplitude  oscillations  in  the  clzunp  ring  when 

mounted  in  a vertical  position. 

The  target  consisted  of  a type  304  stainless  steel  circular 

-2 

plate,  thickness  19.5  mil  (4.95  x 10  cm).  The  density  of  the 

target  material  is  7.788  G cm  total  mass  260.05  G.  This  target 

material  was  chosen  for  two  reasons;  it  is  reasonably  typical  of 

materials  used  in  the  aerospace  industry,  and  the  thermal  diffusivity 

-2  2 -1 

is  sufficiently  low  (4.1  x 10  cm  sec  ) so  that  the  thermoelastic 
effect  is  reasonably  large. 

The  experimental  apparatus  utilized  for  this  Investigation  is 

outlined  in  Fig.  10.  The  laser  used  was  a Holobeam  model 

630-Q  NdjGlass  system.  This  laser  produces  an  output  pulse 

(in  the  Q switched  siode)  of  approximately  3 joules,  with  a pulse 

-9 

width  of  approximately  40  x 10  sec.  The  output  beam  is  about  1cm 
in  diasMter,  with  a beam  divergence  of  2 m rad,  at  an  output  wavelength 
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of  1.06U.  The  beam  is  reflected  from  a Cu  mirror  so  that  it  irradiates 

the  target  at  perpendicular  incidence.  The  direct  beeun  from  the  laser 

2 

has  an  intensity  of  the  order  of  100  MW/cm  , which  is  more  than 
sufficient  to  vaporize  the  target  surface,  which  was  not  desired  in 
this  experiment.  Consequently,  the  beam  was  expanded  using  a negative 
focal  length  lens  to  the  point  below  which  surface  vaporization  occurs. 
In  order  to  insure  that  the  surface  irradiation  was  axially  symmetric, 
the  beam  was  directed  after  expansion  through  a circular  aperture 
which  defines  the  area  of  the  plate  irradiated  by  the  laser,  and 
insured  that  the  laser  irradiated  the  target  at  its  center.  This  was 
done  in  order  to  inhibit  the  generation  of  non-axially  syinnetric  mode 
oscillations  in  the  target. 

The  beam  splitter  in  the  beam  served  to  deflect  a small  portion 
of  the  input  beam  into  a Korad  Photodiode  detector  so  that  the  laser 
pulse  energy  and  power  can  be  monitored.  The  laser  Q switched  pulse 
can,  of  course,  be  varied  in  energy  by  varying  the  capacitor  discharge 
voltage,  but  this  is  accoeq>anied  by  a change  in  the  output  pulse 
width.  Consequently,  the  experiments  were  performed  at  a constant 
discharge  voltage  corresponding  to  an  output  energy  of  approximately 
2j.  Under  these  conditions  the  beam  power  density  %fas  held  to  a 
point  such  that  no  surface  plasma  ms  generated  nor  melting  induced, 
but  such  that  the  surface  teiqperature , based  upon  Bechtel's  treat- 
ment [3],  ms  calculated  to  be  close  to  the  melting  point. 

It  ms  considered  desirable  to  investigate  the  plate  dynamic 


response  resulting  from  the  thermoelastic  stress  utilizing  a method 
which  did  not  require  physical  contact  with  the  plate,  since  such 
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contacts  introduce  boundary  conditions  which  may  perturb  the  system 
dynamic  response.  Since  the  absolute  2ui^lltudes  of  the  stress  waves 
generated  by  the  laser  pulse  are  difficult  to  predict,  owing  to  un- 
certainties in  such  parameters  as  the  reflection  coefficient, 
specific  heat,  and  thermal  conductivity  as  a function  of  surface 
temperature,  it  was  desirable  to  provide  a detection  system  capable 
of  detecting  very  small  amplitudes  of  the  resulting  plate  oscillations. 

Consequently,  the  plate  motion  induced  by  the  laser  pulse  was 
initially  detected  using  an  optical  interferometer.  However,  since 
the  amplitude  of  the  induced  motion  was  in  fact  relatively  large,  a 
much  simpler  detection  method,  as  shoim  in  Fig.  10  was  used.  This 
consists  of  a Bruel  6 Kjaer  model  NM004  capacitive  transducer  and 
associated  circuitry.  The  amplified  signal  is  displayed  on  a 
Tektronix  545B  oscilloscope  and  photographed.  The  detector  was 
mounted  at  a separation  of  lam  from  the  rear  surface  (away  from  the 
laser),  at  the  center  of  the  plate. 

The  oscilloscope  is  triggered  from  a signal  derived  from  the  Q- 
switch  Fockels  cell  trigger,  either  in  the  normal  sweep  or  delayed 
s%ieep  mode.  In  this  manner,  the  plate  motion  is  known  in  time  with 
respect  to  the  occurrance  of  the  laser  pulse,  and  may  be  examined  at 
any  time  during  the  oscillation  utilizing  the  variable  time  delay  of 
the  delaying  sweep. 

Knowing  the  detector  sensitivity,  the  plate  surface  to  detector 
•eparation,  the  asqplifier  gain,  the  capacitor  polarization  voltage, 
and  the  oscilloscope  vertical  sensitivity,  the  overall  sensitivity  of 
the  detector  is  known.  Therefore,  the  amplitude  of  the  plate 
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oscillation  can  be  directly  determined  from  the  oscilloscope  traces. 

In  our  case,  the  total  sensitivity  of  the  detection  system  is  1.905  x 
10  ^ volts/ln.  The  maximum  anplltude  of  the  Laser  induced  thermo- 
elastic  oscillations  observed  in  this  experiment  was  1.3  x 10  ^ in, 
or  roughly  6.7%  of  the  plate  thickness.  Under  these  conditions,  one 
is  confident  that  the  general  assumptions  of  plate  theory  are  valid. 

The  temperature  rise  at  the  rear  surface  of  the  target  plate 
was  monitored  using  a thermistor  detector.  Since,  for  the  values 
of  plate  thickness,  thermal  diffusivity,  and  pulse  width  used  in 
this  experiment  the  plate  is  thermally  thin,  the  temperature  at  the 
rear  surface  is  a valid  mteasure  of  the  total  energy  ^d>sorbed.  The 
data  then  indicate  an  absorbed  energy  of  1.16  joule,  which  is 
consistent  with  estimates  of  the  reflection  coefficient  at  1.06m  and 
the  reflection  coefficient  of  the  Cu  miirror.  This  value  of  absorbed 
energy  then  yields  an  estimate  of  the  surface  temperature  of  1 x 10^  °C. 
which  is  significant  only  in  that  it  agrees  with  the  non-vaporization 
or  melting  assumption.  It  is  not  used  in  our  analysis  of  the  data, 
since  we  are  focusing  our  interest  on  the  plate  oscillation  induced 
by  the  laser  Impulse. 

The  general  features  of  the  plate  dynamic  response  are  illus- 
trate by  the  oscilloscope  traces  shown  in  Figs.  11-14.  These 
figures  are  oscilloscope  traces  for  various  sweep  speeds  of  the  plate 
oscillations  induced  by  the  laser  pulse.  Several  general  features 
of  the  plate  response  can  be  imamediately  seen  in  these  traces.  First, 
the  initial  response  of  the  plate  consists  of  a rapid  mntion  in  a 
direction  toward  the  face  of  the  plate  Irradiated  by  the  laser. 
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FIGURE  12 

Sweep  Speed  10  msec/cm. 


FIGURE  13 

Sweep  Speed  5 msec/cm 


FIGURE  14 

Sweep  Speed  2 msec/an 


This  would  be  expected  if  the  primary  excitation  of  the  plate  were 
by  thermoelastic  effects.  An  impulse  arising  from  surface  blow  off, 
plasma  formation,  or  thermal  effects  arising  from  expansion  of  the 
air  in  contact  with  the  heat^ad  surface  would  be  expected  to  induce 
initial  motion  in  the  opposite  direction.  Second,  the  induced  plate 
oscillations  contain  many  higher  order  modes  which  have  been  included 
in  the  theoretical  model  of  this  process.  It  is  clfar  that,  as 
expected,  these  higher  order  modal  frequencies  are  dampel  at  varying 
damping  rates,  so  that  after  a time  of  the  order  of  120  m sec,  the 
plate  can  be  seen  to  be  oscillating  in  its  fundeunental  (l.e.,  lowest 
order  mode)  frequency.  As  mentioned  previously,  we  have  not  modeled 
these  variable  damping  rttes  in  the  theoretical  treatment  of  this 
problem.  We  have  made  use  of  the  fact  that  the  plate  ultimately 
oscillates  in  its  lowest  order  mode  to  obtain  the  necessary  material 
constants  for  use  in  the  theoretical  model,  as  discussed  in  Appendix  II. 

Third,  the  general  response  of  the  plate  consists  of  the  initial 
tremsient,  the  fundamental  ^md  higher  order  mode  oscillations,  and, 
in  addition,  a superlo^sed  exponential  relaxation  response,  asso- 
ciated with  an  effective  time  constant  of  the  order  of  17  ■ sec. 

That  this  portion  of  the  plate  dynamic  response  is  associated  with 
laser  heating  is  shown  by  the  fact  that  it  is  absent  vdien  the  plate 
oscillations  are  Induced  by  mechanical  impulse.  This  is  discussed 
in  Appendix  II.  It  is  attractive  to  associate  this  effect  with  a 
possible  theraal  pre-stress  condition,  and  experiments  relating  to 
that  effect  are  continuing.  To  date,  however,  vre  have  not  Included 


a Bodsl  of  this  affect  in  our  theoretical  treatawnt 
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The  comparison,  then,  between  the  theoretical  calculations  of 
the  plate  dynamic  response  emd  the  experimental  results  is  performed 
in  the  following  manner.  The  theoretical  response  curve  (Fig.  5)  is 
based  upon  material  constants  chosen  so  as  to  agree  with  the  actual 
observed  fundamental  mode  frequency,  as  discussed  in  the  Appendix  II. 

On  this  curve,  the  point  on  the  time  scale  (1.0)  then  corresponds  to 
a time  of  8.33  m sec.,  and  the  other  time  points  are  correspondingly 
known.  The  theoretical  curve  (Fig.  5)  is  then  mounted  on  a flat 
surface,  and  a transparency  made  from  the  oscilloscope  trace  is 
projected  onto  the  theoretical  curve.  The  magnification  of  the 
optical  system  is  then  varied  so  that  the  time  scale  of  the  oscillo- 
scope trace  agrees  with  the  horizontal  (time)  scale  of  the  theoretical 
curve.  Then,  maintaining  the  time  scale  fixed,  the  flat  surface  on 
which  the  theoretical  curve  is  plotted  is  simply  rotated  so  that  the 
vertical  (amplitude)  scales  are  normalized  at  any  arbitrary  point. 

The  projection  of  the  experimental  curve  is  then  simply  plotted 
directly  onto  the  theoretical  curve,  yielding  the  coaq>arison  between 
theory  and  experiment  shown  in  Fig.  15.  This  figure  indicates  the 
comparison  in  response  over  a time  corresponding  to  two  fundamental 
periods,  or  about  17  m sec.  Comparison  of  the  theoretical  curve  and 
the  experimental  data  taken  at  a delay  time  of  7 m sec,  using  the 
variable  delay  of  the  triggered  sweep,  is  sho%m  in  Fig.  16. 

Since,  as  discussed  previously,  we  have  not  included  a consideration 
of  the  variable  higher  order  mode  damping  in  our  model,  it  is  not 
particularly  useful  to  extend  our  comparisons  much  beyond  a time  scale 
of  20  a sec.  In  fact,  it  is  clear  from  the  coaq;>arlson  of  Fig.  15  that 


lA 


Vertical  Scale  - Experimental  Data  expanded  by  x 2. 
Theoretical  Scale  Arbitrary 

Theoretical  Data 

Experimental  Data 
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several  of  the  higher  order  mode  frequencies  have  already  sufficiently 
decayed  and  they  are  not  contributing  significantly  to  the  actual 
plate  response. 

Using  the  comparisons  illustrated  in  Fig.  15  and  16,  it  seems 
clear  that  the  general  features  of  the  response  are  predicted  quite 
reliably.  We  believe  that  the  experimental  data  can  be  taken  as  a 
verification  of  the  thermoelastic  model  of  the  laser  interaction  at 
this  power  level. 
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8.0  CONCIiJSIONS 

(a)  A mathematical  model  of  the  Laser  Induced  thermoelastic  motion 
of  a plate  was  derived.  The  model  consists  of  three  partial 
differential  equations  with  associated  boundary  conditions. 

Two  of  these  equations  decouple  and  characterize  plate  flexure. 
The  remaining  equation  describes  plate  stretching. 

(b)  Laser  Induced  heating  is  Incorporated  Into  the  model  by  using 
Bechtel's  surface  heat  generation  model  (3]. 

(c)  The  equations  were  used  to  solve  the  problem  of  a thin  circular 
plate  clamped  at  its  boundary.  The  circular  plate  is  subjected 
to  Laser  Irradiation  at  its  center,  and  normal  to  the  plate 
surface. 

(d)  Laboratory  experiments  were  performed  to  measure  the  center 
deflection  of  the  plate  under  Laser  irradiation.  The  results 
of  these  experiments  conqpare  well  with  those  obtained  by 
calculations  using  the  mathematical  model. 
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APPENDIX  I:  NUMERICAL  EVALUATION  OF  g. (t) 
1 

In  order  to  resolve  the  forced  transverse  motion  problem  In 
section  6.1  it  is  necessary  to  evaluate  the  function 

91  (ti)  - M(tj^)-2nR^  M(e.)sin  27rR^  (tj^-0)d0  (A-1) 

o 

where  M(t^)  is  defined  by  (93) . To  simplify  the  forthcoming  analysis 
%#e  will  omit  the  subscript  i on  g^,  1 on  t^^  and  let  (U“2tiR^.  Then  (A-1) 
can  be  re%n:itten  as  follows: 
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where , 


2)c 

AF  - u)  / M(e)co8  ojeae 

2(k-l) 


2k 

AG  = u)  / M(e) 

^ t 

2(k-l) 


sin  (D6d6 


(A-9) 


(A-10) 


and  P -G  =0. 
o o 


To  evaluate  the  integrals  (A-9)  and  (A-10) , we  use  a second  order 
Lagrange  interpolation  formula  to  approximate  M(6)  on  the  interval 

^^2(k-l)'  *^2k^' 


M(0)  - y L +Y  L +YL 
' 2k-2  2k-2  2k-l  2k-l  2k  2k 


where 


<A-11) 


(A-12) 


Sk  - <«-'=2k-2"»-'2k-l>''='''''  > 


where,  Y^“M(t^)  j-0,1, . . . ,2N. 


Substitution  of  (A-11)  and  (A-12)  into  (A-9)  and  (A-10)  yields 
the  results 


* 2il  ^\^2k'*^\^2k-2^'^  „2  ^®2k-2”®2k^ 


k-1  2k  2k  2k-2  2k-2 


(A-13) 


^k-1  “ ’^2k-2^2k-2"’^2k^2k  * 2S  ^*k®2k'*‘**k®2k-2^  * ^2  ^Sk"^2k-2^ 


(A-14) 
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where , 


fi  = ujAt 


3Y  -4Y  +Y 

2k-l  2k-2 

Y -4Y  +3Y 

2k  2k-l  2k-2 

Y -2Y  +Y 

2k  2k-l  2k-2 


(A-15) 


S . = sin  a)t„,  = sin  2kfi 
2k  2k 

C^,  * cos  a)t_,  = cos  2kfi 
2k  2k 


Consequently, 

= Y 


2k'  *2k  V2k‘*^kSk  ' N. 


(A-16) 


To  check  the  accuracy  of  this  approximation,  the  following  example, 
having  a closed  form  solution,  was  tested. 


M(t)  = e 
In  this  case 


-at 


g(t) 


2 -at 

g e -<ti(a  sin  tot -tocos  a)t) 

2 2 
a + u) 


(A-17) 


(A-18) 


Table  (A-1)  shows  the  comparison  of  the  exact  result  from  (A-18)  with 

the  approximation  (A-16)  for  the  case  t •!,  M«10,  a^27T,  a>1.84. 

m&x 

The  agreement  with  N-10  is  seen  to  be  quite  good.  Further  tests  with 

values  of  U)  up  to  200ir  and  N>10  indicate  that  the  accuracy  does  not 

decrease  significantly  as  a>  increases.  Another  test  with  t ■!, 

n&x 

h)*2irR^^,  where  N>100  was  run  and  the  agreement  betmen 

the  exact  and  approximate  values  was  at  least  seven  aitd  sometimes 
eight  decimal  places. 


Based  on  these  observations,  the  final  calculations  in  section 


(6.1)  were  run  using  t -2  and  N«'200. 

nax 


APPENDIX  II;  EXPERIMENTAL  DETERMINATION  OF  PLATE  STIFFNESS. 


It  became  clear  even  during  the  Initial  stages  of  this  experi- 
ment, in  which  interferometric  techniques  were  utilized  to  study  the 
plate  dynamic  motion,  that  the  fund2unental  frequencies  observed  were 
in  general  significantly  different  from  those  predicted  by  h2mdbook 
values  of  the  material  constemts  for  our  selected  target  material. 

It  was  therefore  considered  necessary  to  perform  some  calibration 
experiments  relating  to  the  determination  of  the  effective  plate 
stiffness,  which  determine  the  predicted  plate  response  to  the  laser 
pulse.  Several  related  experimental  techniques  were  used  in  this 
regard,  namely  excitation  of  the  plate  oscillations  by  mechanical 
impulse,  electromagnetic  excitation  of  plate  oscillations,  holo- 
graphic analysis  of  static  deflection,  and  direct  measurement  of 
the  purely  fundamental  mode  frequencies  during  the  laser  excitation. 

Since  this  latter  method  is  most  directly  related  to  the  main 
experimental  approach,  we  will  describe  it  first.  The  circuit 
outlined  in  Fig.  A1  is  used.  The  signal  from  the  photodiode 
detector  triggers  the  delaying  sweep  of  the  545B.  After  a time 
determined  by  the  variable  delay,  the  main  sweep  of  the  545B  is 
triggered,  and  the  main  sweep  signal  serves  as  enable  gate  to  the 
IC  gate  circuit.  When  this  gate  signal  is  present,  the  detected 
signal  from  the  capacitive  transducer  is  passed  to  the  5245  electronic 
counter  operated  in  the  period  mode.  Each  time  the  laser  is  fired,  a 
measurement  Is  obtained  of  the  fundamental  plate  frequency  by 
delaying  the  main  sweep  trigger  until  the  higher  order  modes  have 
damped  out.  The  period  observed  in  this  manner  was  then 
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(8.2630  ± 0.0277)  x 10  ^ sec.,  corresponding  to  a fundamental  frequency 
of  (121.02  ± 0.58)  Hz. 

The  response  of  the  plate  to  a mechanical  impulse  was  investigated 
using  the  apparatus  outlined  in  Fig.  A2.  In  this  experiment,  the 
mechanical  Inqpulse  was  produced  by  allowing  a small  copper  sphere  to 
strike  the  plate  at  its  center.  Since  the  target  plate  is  mounted  in 
a horizontal  position,  as  discussed  in  section  7.0  of  this  report, 
the  sphere  is  suspended  by  elastic  supports  so  that  it  does  not 
contact  the  plate  in  its  equilibrium  position.  A signal  derived  from 
the  interruption  of  a He-He  laser  beam  passing  parallel  to  the  plate 
face  was  used  to  trigger  the  delaying  sweep  of  the  oscilloscope  in 
the  scune  manner  as  the  apparatus  used  in  the  main  experiment. 
Oscillogrcuns  of  the  plate  response  to  this  excitation  are  illustrated 
in  Figs.  A-3  through  A-5. 

Several  conclusions  can  be  reached  from  these  responses.  First, 
it  is  clear  that  the  Initial  plate  response  is  in  a direction  away 
from  the  mechanical  impulse,  as  would  be  expected.  This  is  in  fact 
the  opposite  effect  from  the  thermoelastic  response  to  the  laser 
inpulse.  Second,  the  oscillations  induced  by  the  mechanical  inpulse 
are  symnetric  about  the  baseline.  This  is  in  contrast  to  the  laser 
induced  plate  response  which  exhibits  am  exponential  transient  in 
addition  to  the  plate  oscillation.  We  therefore  believe  that  this 
transient  effect  is  a result  of  the  thermoelastic  excitation,  and 
is  not  characteristic  of  the  plate  itself.  In  addition,  the  frequency 
spectrum  of  the  response  to  this  method  of  excitation  contains  fewer 
higher  order  modes  than  is  the  case  for  thermoelastic  excitation. 


Elastic 

SupportX  /Metal  Ball 
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Figure  A2:  Mechanical  Impulse  Excitation 
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FIGURE  A-3 

Plate  Response  to  Mechanical  Impulse 
10  ms/div  5v/div 
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FIGURE  A-4 
5 ras/div  5v/div 


f 


FIGURE  A-5 
2 ms/div  5v/div 
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It  is  clear  that  they  are  damped  at  variable  rates,  as  observed  in  the 
laser  excitation.  Lastly,  it  can  be  seen  that  the  fund^ullental  frequency 
is  in  agreement  with  the  laser  excited  frequency,  as  one  would  expect. 

In  order  to  study  the  plate  vibration  frequency  spectrum,  the 
apparatus  outlined  in  Fig.  A6  has  been  used.  In  this  experiment, 
the  plate  oscillations  aure  excited  by  means  of  an  electromagnetic 


transducer  (B&K  Model  0002)  driven  by  a varied>le  frequency  oscillator. 
The  voltage  derived  from  a 390  n resistor  in  series  with  the  transducer 
is  used  in  the  feedback  mode  to  maintain  a constant  amplitude  of  the 
plate  excitation  as  the  frecjuency  is  varied.  This  is  necessary  since 


the  transducer  presents  an  Inductive  reactance  load  to  the  oscillator. 
Using  this  apparatus,  the  various  modal  oscillation  frequencies  may 
be  separately  excited  by  the  oscillator.  The  first  five  modes  were 
easily  observable  in  this  test.  The  observed  frequencies,  at  a plate 
teii$>erature  of  76”F  were  120.1,  467,  1055,  1858,  and  2900  Hz. 

The  value  of  Young's  modulus  for  the  material,  or  perhaps  more 
properly,  the  effective  plate  stiffness,  may  be  determined  from  these 


data  using  the  relation: 


1 1 

i-na^  Li2(1-v^)i;  J 


(A-19) 


where 


- modal  eigen  value 
h > plate  thickness 

a " plate  (free)  radius 

E « Young's  Modulus 
C ■ plate  density 


V > Poisson's  ratio 


Figure  A6:  Electromagnetic  Excitation  of  Plote 


I 


I 

I 86 

r 

During  the  course  of  this  experiment,  it  became  clear  that  the 
resonant  frequency  observed  was  in  fact  highly  sensitive  to  the  ambient 
plate  ten^erature.  In  general,  it  has  been  observed  that  the  lowest 
order  mode  frequency  shifts  at  a rate  of  4.88  Hz/'F.  Consequently, 
the  data  have  been  taken  at  or  near  the  76°F  value  mentioned  above. 

We  presently  feel  that  this  effect  is  probably  due  to  a thermal 
pre-stress  condition  eurising  from  the  fact  that  the  thermal  masses 
of  the  massive  clasp  and  thin  plate  are  greatly  different.  Therefore, 
ambient  temperature  chamge  affects  the  plate  more  rapidly  than  the 
clamp.  Theoretical  analysis  of  this  effect  is  continuing  at  present, 
along  with  further  experimental  investigation  of  the  effect  of  the 
tenperature  difference  (between  clamp  and  plate)  on  the  resonant 
frequencies. 

In  addition,  it  was  considered  desirable  to  investigate  the 
static  deflection  characteristics  of  the  plate  as  em  independent 
study  of  the  material  constants  of  the  plate.  The  apparatus  used  is 
sho%m  in  Fig.  A7.  Weights  of  10,  20,  50,  and  100  G were  applied  to 
the  plate,  and  double  exposure  holograms  %rere  made  of  the  plate. 

That  is,  a hologram  of  the  undeflected  plate  (no  load)  was  made  by 
exposing  the  holographic  plate,  the  holographic  plate  was  covered, 
the  load  applied,  and  a second  exposure  was  made.  After  the  hologram 
was  developed  and  dj led,  it  was  replaced  in  the  plate  holder  and 
illuadnated  lay  the  reference  beam  only.  A photograph  was  then  made 
of  the  holographic  image  thus  formed.  This  photograph  then  shows  an 
image  of  the  plate,  superisiposed  on  vdiich  is  a series  of  interference 
fringes  resulting  from  the  static  deflection  of  the  plate.  Displacement 
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values  are  then  determined  by  counting  the  interference  fringes  from 

% 

the  undeflected  edge  of  the  plate  image  to  emy  point  on  the  plate. 
Since  each  fringe  shift  corresponds  to  a phase  shift  of  X/2,  where  X 
is  the  laser  wavelength,  and  therefore  corresponds  to  a total  plate 
deflection  of  X/4,  the  displacement  ceui  be  determined  to  a precision 
of  X/4,  approximately,  after  correction  for  the  incidence  and  reflec- 
tion angles.  Since  the  deflections  are  then  determined,  the  plate 
characteristics,  stiffness  and  Young's  modulus,  may  be  determined  in 
two  separate,  but  related,  methods  of  analysis.  First,  the  equation 
of  the  static  deflection  of  a centrally  loaded  plate  is  given  by: 


e A8:  Static  Deflect! 
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This  procedure  was  carried  out  for  the  target  plate  used  in  the 
laser  irradiation  experiments,  and  for  a non-irradiated  sheet  cut  from 
the  same  original  sheet. 

The  results  indicated  that  the  values  of  E determined  from 

equations  {A-21)  and  (A-22)  agree  to  within  1%.  However,  the  value 

of  E determined  in  this  manner  for  the  non-irradiated  S2ui^le  was 

€*  2 

unreasonably  large  (of  the  order  of  80  x 10  Ib/in  ) . Since  the 
values  obtained  from  the  two  methods  of  analysis  agree,  and  since 
the  ratio  of  the  exp>erimental  value  of  the  slope  in  equation  (A-3)  to 
the  experimental  value  of  the  slope  in  equation  (A-4) , that  is: 


(A-23) 


r-a/e 


agrees  to  within  1.5%  of  the  expected  value,  it  is  believed  that  the 
radial  dependence  as  given  in  equation  (A-23)  is  correct.  It  is  noted 
that  the  experimental  deflections  in  this  test  were  intentionally 
limited  to  small  values  (1.6  x 10  ^ in,  maximim) . It  is  possible 
that  the  large  values  observed  in  this  test  are  an  additional 
indication  of  a thermal  pre-stress  condition  referred  to  previously. 

A more  interesting  observation,  however,  was  that  the  value  of 
E determined  for  the  laser  irradiated  plate  was  about  1.7  times  larger 
than  the  non-irradiated  plate.  During  the  initial  measurements  on 
the  laser  irradiation  results,  the  plate  was  of  course  irradiated  many 
tisMS,  and,  because  the  beam  was  not  initially  expanded,  a laser 
produced  surface  plasma  %«as  Induced  several  times.  Subsequently,  the 
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plate  was  pulsed  repeatedly  at  lower  power  density.  Consequently,  we 
do  not  know  whether  the  apparent  Increase  was  due  to  the  initial  high 
intensity  incises,  the  later  lower  Intensity  impulses  or  both. 
Furthermore,  after  the  second  target  plate  had  been  used  for  the 
thermoelastic  study  euid  the  required  data  obtained,  we  irradiated  it 
at  higher  power  densities,  and  utilized  the  electromagnetic  trauisducer 
(as  outlined  above)  to  determine  the  fundamental  resonant  frequency. 

An  increase  in  the  value  was  Immediately  observed,  again  raising  the 
question  of  whether  the  high  power  pulse  has  resulted  in  an  effective 
increase  in  plate  stiffness.  However,  until  the  effects  of  the 
thermal  pre-stress  condition  are  better  understood,  we  believe  it 
would  be  premature  to  assign  the  observation  to  a change  in  the 
material  constants  of  the  plate. 


